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ABSTRACT 



Simplified finite element modeling of a stiffened cylinder subjected to 
underwater explosion was investigated. The use of smearing the stiffeners into the 
base structure as well as beam modeling using SOR (Surface of Revolution) beam 
elements were used in the model simplification process. The dynamic response and 
overall global deformation were then compared between the fully discretized stiffened 
cylinder model and the simplified models. The study first examined the 
effectiveness of smearing stiffeners onto a flat plate. The smearing of stiffeners into 
a cylindrical shell orthotropically was then examined. Next, beam modeling of both 
unstiffened and stiffened cylinders was investigated. Finally an integrated beam/shell 
model of a stiffened cylinder was developed. These models were subjected to the 
same underwater explosive loading for numerical study. The analysis showed that 
when comparing the dynamic responses caused by underwater explosions between 
the discrete model, the beam model, and the beam/shell model of a stiffened cylinder, 
globally similar results could be produced. 
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I. INTRODUCTION 



Detailed finite element modeling of a submerged structure including internal 
structures subjected to underwater explosion is costly. The cost can be equated in both 
computational time and more importantly manpower. As a result, a simple, but still 
accurate model of a structure was investigated. The objective of the study was to develop 
an efficient means of representing the stiffened cylinder while still providing accurate 
dynamic response predictions to underwater explosion. Past work at the Naval Postgraduate 
School focused on comparison of numerical modeling to the physical testing. Cylindrical 
models have been previously subjected to various types of underwater explosions and their 
results have been closely approximated by the computer code. Fox, Kwon, and Shin 
[Ref. 1] and Chisum [Ref. 2] have all demonstrated that the finite element and boundary 
element programs closely approximate experimental results. Their research work clearly 
indicates that limited numerical modeling can be conducted without the need to physically 
build the actual structure for experimental testing. 

The studies [Ref. 1] did not explore modeling cylindrical shells with internal 
structures. This type of modeling involves adding internal equipment and structures to the 
outer shell, which dramatically increases the degrees of freedom of the problem. This 
directly correlates to a long computational time for solutions. This study will examine 
two means by which a simplified model and discrete model can be numerically tested to 
produce similar dynamic responses due to an underwater explosion. The first technique 
models the stiffened cylinder with SOR(Surface of Revolution) beam elements and the 
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second technique involves smearing the stiffeners onto the cylindrical shell structure. 

Using a material smearing technique is prevalent in composites where smearing 
material properties reduces the complexity of the structure when subjected to various forms 
of loading. Raftenberg [Ref. 3] demonstrated that smearing in structures required that the 
structure be entirely axiymetric or have regions which varied circumferentially about its 
axis of symmetry with rapid periodicity. The smearing technique was used for static loading 
only. Pitarresi, Caletka, Caldwell , and Smith used a similar finite element material-smearing 
technique to dynamically test a printed circuit board with favorable results [Ref. 4]. Both 
studies dramatically reduced both complexity and computational time with their 
simplification methods. 

The finite element and boundary element programs utilized for providing the 
numerical analysis for this study are called USA/DYNA3D. VEC/DYNA3D is a finite 
element analysis code developed by Lawrence Livermore National Laboratory [Ref. 5] . The 
Underwater Shock Analysis, USA [Ref. 6 ] is a boundary element code that is based on the 
Doubly Asymptotic Approximation (DAA) which was developed by Geers [Ref. 7]. The 
linkage between the two codes was developed in 1991 and has provided acceptable results 
of dynamic responses of cylindrical shells. 

One of the applications of the USA code developed by DeRuntz [Ref. 6] allows the 
beam to be modeled with SOR (Surface of Revolution) elements. The SOR element 
approximates a beam element with a surface of revolution element. The SOR element 
contains a number of sub-elements around the circumference of each surface of revolution 
element. The SOR branch interpolates sub-elements using a circumferential harmonic 
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function which ultimately reduces the number of surface of revolution degrees of freedom. 

The original discrete model of the cylinder analyzed in this study was provided by 
Protection & Weapons Effects Department, Code 67 Carderock Division of the Naval 
Surface Warfare Center and is provided in Figure 1. This cylinder is called the ONR 
cylinder and is a test vehicle which is currently being planned for a series of underwater 
explosions. This study will first examine the effects of smearing stiffeners onto a flat plate. 
Next, the discrete model will be modeled without internal stiffeners with Belytschko-Lin- 
Tsay 4 node shell elements due to its computational efficiency in DYNA3D program. The 
unstiffened model will be modeled as a beam with SOR elements and an integrated 
beam(SOR)/shell model will be developed. All three models will be subjected to simulated 
underwater explosions and their responses will be analyzed and discussed. 

The discrete model with internal stiffeners will then be modeled. Based on the study 
of smearing stiffeners into a flat plate, the ONR stiffeners are smeared into the shell of the 
cylinder orthotropically. Beam model with SOR elements will be developed as well as 
integrated beam/shell model. The internal stiffeners will be incorporated into beam model 
by adjusting section properties of beam. The discrete and two simplified models will be 
numerically analyzed by applying the explosive charge first at the symmetric point of all 
three models. Detailed analysis will be provided at beam/shell interface as well as 
endplates and deep frame stiffeners to compare discrete model to simplified models. The 
three models will then be exposed to an off-center charge and analysis will be performed as 
previously noted. 



3 




Figure 1. Naval Warfare Surface Center's ONR Cylinder 
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n. UNDERWATER EXPLOSION CHARACTERISTICS 



The underwater explosion testing that is simulated by USA computer code provides 
the user many options. For this study, however, the same options and explosion 
characteristics were used so that a simlilar comparison between the discrete or full model 
and simplified models could be made. All models developed in this study were subjected 
to an underwater shock wave modeled as a spherical wave with a peak amplitude of 2585 
psi and a decay constant of 0.3493 msec. This approximation to the actual shock wave uses 
an exponential pressure relationship and is good for pressure greater than about 1/3 of the 
peak value [Ref. 8]. This peak pressure pulse is generated using 60 lbs charge of HBX-1 
explosive at a 10 foot depth and a standoff distance of 27.2 feet. The attack geometry is 
provided in Figure 2. The pressure profile is provided in Figure 3 for a reference. 

The 10 foot underwater depth was selected in order to fall within the gas bubble 
radius of 15.8 feet. This allows the gas bubble to vent to the surface and not migrate to the 
surface of the cylindrical surface. If the bubble were not to vent to the surface then a 
significant amount of explosive energy would be subjected to the structure(approximately 
47%) [Ref. 8]. This was done in order to eliminate the effects of bubble pulse on the 
simulations perfomed in this study. Free surface effects were also ignored in this study. 
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Figure 2. Attack Geometry 



60# HBX-1 UNDEX Pressure Profile 




Figure 3. Pressure Profile for HBX-1 60 Pound Charge 
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ID. FLAT PLATE SMEARING 



The first area of model simplification that was studied involved examining the effects 
of smearing stiffeners onto a flat plate. The aim of the study was to find an optimal way to 
smear the stiffeners without affecting the global dynamic response of the flat plate. From 
this study, the objective was to use these results for finding an efficient means of smearing 
the stiffeners onto the cylindrical shell. 

The plates used in this study were either clamped or simply supported at the edges 
and were subjected to a unit step pressure load of lOOOpsi. All plates examined were 88" 
x 88" with a stiffener spacing of 4 inches. The stiffeners are 1.25" thick and 0.25" wide and 
were modeled in both directions. A comparison between the discrete and smeared models 
was accomplished using VEC/DYNA3D finite element programs. The discrete model had 
the stiffeners fully modeled to the base plate whereas the smeared model had the stiffeners 
smeared onto the base plate. The smeared value was simply an adjustment in the material 
property of the base material. The adjustment is mad within the material definition portion 
of INGRID. The smeared value was obtained using a well known formula derived from 
Timoshenko and Woinowsky-Kreiger [Ref. 9]: 



E' = E 



12 El t (1-v 2 ) 

ah 3 
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where 



E'= fictitious smeared elastic modulus 
E = elastic modulus of plate and stiffener 
f. = moment of inertia of the stiffener 
a = spacing between two stiffeners 
h = thickness of the plate 
v = Poison's ratio 

The density of the base material was also adjusted to take into account the mass of the 
stiffeners. 

The discrete model of a 0. 1875" thick plate with clamped edges was generated and 
subjected to the unit step pressure loading. Figure 4 is the dynamic response of the discrete 
model at 0.5 msec. The transverse deformation plot is magnified by a scale factor of 10 for 
clearer visualization of the deformation. Next the smeared value was computed for a 0. 1875 
" thick plate with the given parameters for the stiffeners. The smeared model was subjected 
to the same pressure loading and Figure 5 shows the deformation of the smeared model. The 
fictitious smearing value was approximately 85 times larger than the original material's 
elastic modulus ( E's 85*E). 

< 

This modeling shows the case where the stiffeners dominated the overall response 
of the original structure. Notice that when comparing the two results, it was observed that 
the discrete model had a much larger deformation gradient near edges of the plate. The area 
near the edges of the discrete model governs the global deformation of the structure and for 
discussion purposes will be called the critically deformed area. The glaring difference in 
global deformation between the discrete and smeared models can be accounted for by the low 
stiffness of the discrete model base plate. In other words, the original discrete model has 
1/85 of the globally smeared model's stiffness in the base plate. As a result, stiffeners in the 
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Figure 4. Global Deformation Plot of Discrete Model of 0.1875" Plate With Clamped 
Edges, Magnified by a factor of 10, taken at 0.5msec 




Figure 5. Globally Smeared Model of 0.1875" Plate With Clamped Edges 
Magnified by a factor of 10, taken at 0.5msec 
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mid-section of the discrete model's base plate are less constrained from the side edge 
constraints. This produces a very large gradient near the edges of the base plate and a 
minimum gradient at the mid-section of the plate. The contrapositive of this analysis is that 
the smeared plate is uniformly constrained in both directions therefore the globally smeared 
plate has a uniform deformation gradient over the entire plate. Although in Figures 4,5 the 
deformation is shown at 0.5 msec, the comparison is very similar even at 3.0 msec. 

The most effective way of simplifying the discrete model and obtaining the proper 
deformation gradient at the critically deformed area is to smear locally. That is, leave the 
discrete stiffeners modeled around the edges of the plate in the critically deformed area, 
however locally smear the remaining stiffeners into the base plate as in the globally smeared 
model using the same smearing value as in the globally smeared model. Figure 6 shows the 
dynamic response of the locally smeared versus the discrete model. The local smearing did 
not affect the critically deformed area. As a result the dynamic response of the locally 
smeared plate was in excellent agreement with the discrete model. 

A second study of the same stiffened plate was undertaken using simply supported 
boundary conditions at the plate edges. The dynamic responses produced similar results to 
what was observed for the clamped plate modeling. The simply supported 0.1875" thick 
plate models (discrete and globally smeared) are provided in Figure 7. The results were 
very similar out to 3.0 msec so only 0.5 msec is shown. 

Next, a study was conducted to analyze the case where the ratio of fictitious elastic 
modulus to original elastic modulus is lowered. By increasing the plate thickness to 1 ", the 
ratio was reduced to 2 (E® 2E). This time, the results were expected to show excellent 
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6a) Locally Smeared Model 




6b) Discrete Model 

Figure 6a, b. Comparison of Locally Smeared model of 0.1875” Plate versus Discrete 
Model of 0.1875” Plate. Both With Clamped Edges and with 
magnification factor of 10, taken at 0.5 msec 
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(a) Discrete Model with Simply Supported Edges 




(b) Globally Smeared Model 

Figure 7. 0.1875" Plates With Simply Supported Edges: (a) Discrete Model (b) 

Globally Smeared Model. Both with 10X magnification , at 0.5 msec 
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agreement between the discrete and globally smeared models due to the increased stiffness 
in the base plate. Again the discrete model was generated and subjected to the same unit 
step pressure loading. Figure 8 shows the result of the 1" discrete plate model versus the 
globally smeared model with both models having clamped edges. Figure 9 compares both 
discrete and globally smeared 1" plates with simply supported edges. As predicted the 
responses were very close to each other for both sets of boundary conditions. 

Based on the study conducted on smearing flat plates, the smearing technique 
directly supported the well known and established formula derived in Ref. 9. The present 
study provided an additional criteria for determining the proper approach to smearing 
stiffeners onto a base structure based on the original plate's thickness and elastic modulus. 
If the ratio between fictitious elastic modulus and original elastic modulus is too high 
erroneous solutions will be obtained if global smearing is undertaken. If the ratio is small, 
then global smearing is an effective means for simplifying the structure. 
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(a) 1” Discrete Model With Clamped Edges 




(b) Globally Smeared 1" Plate Model 

Figure 8. 1" Plates With Clamped Edges: (a) Discrete Model (b) Globally Smeared 
Model. Both with 10X magnification factor and taken at 0.5msec. 
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(a) Discrete Model With Simply Supported Edges 




(b) Globally Smeared Model 

Figure 9. 1" Plates With simply Supported Edges: (a) Discrete Model (b)GIobally 
Smeared Model. Both at 10X magnification factor and at time 0.5 msec 



15 









16 



IV. MODELING OF UNSTIFFENED ONR CYLINDER 



When first considering the task of simplifying a stiffened cylinder, a detailed study 
was undertaken to investigate simplification techniques in general. Two techniques were 
investigated for modeling internal structures such as stiffeners into the external shell. The 
technique previously investigated studied smearing the stiffeners onto the base structure. 
This knowledge will be used when investigating orthotropically smearing the stiffeners into 
the cylindrical shell later in this report. The second simplification technique involves using 
SOR elements to effectively model the cylindrical shell and internal structures. This section 
of the work details how to model an unstiffened cylinder with SOR elements. The SOR 
beam model will be developed and compared to the unstiffened shell. Finally, an integrated 
SOR beam/shell model of the unstiffened cylinder will be generated and all three models will 
be subjected to the same explosive loading. 

A. DISCRETE UNSTIFFENED ONR CYLINDER MODEL 

The three-dimensional model is shown in Figure 10. Its dimensions conform to 
those depicted on Figure 1 except for no internal stiffeners. This particular model uses HY- 
100 steel for the shell and mild steel for the endplates. Both materials for modeling 
purposes were considered an isotropic elastic material. Material properties used for this 
model are provided in Table 1. 
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Mild Steel 



HY-100 



Density 


(lbm/ft 3 ) 


Poisson's Ratio 


Yield Strength 


(ksi) 


Young's Modulus 


(psi) 


Tangential Hardening 


Modulus 


(psi) 



490.0 


490.0 


0.3 


0.3 


32.0 


108.0 


2.9xl0 7 


2.9xl0 7 


5.1xl0 4 


5.02x10' 



Table 1. Material Properties for ONR Cylinder 




Figure 10. ONR Cylinder Modeled Without Internal Stiffeners 
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Although the model shown is a full model, a half model of the cylinder was used in 
the numerical analysis to take advantage of symmetry condition. The yz-plane was the 
symmetric plane from Figure 10. Even with a half cylinder model, there are 20 elements 
around the half circumference and 45 elements in the axial direction with an element length 
of 5.625 inches. The element length represents the spacing between the shallow frame 
stiffeners. 

The discrete model was subjected to an underwater shock with USA/D YNA3D 
programs with shot geometry shown on Figure 2. The major response of the cylinder to the 
shock is in the yz-plane. Since there was a boundary condition imposed on yz-plane, there 
is no global motion in the x-axis. As a result, only global deformation responses in the z- 
direction are compared in Figure 1 1 at discrete states of time. The plots are magnified by 
a factor of 50 for improved visual ease. At early time, the deformation appears symmetric 
except for slight pinching near endplates. The endplates behave similarly to the clamped 
edges on the flat plates in that a large deformation gradient forms around each endplate. 

Another way to analyze the dynamic response is to plot the displacement, velocity 
and acceleration at desired points on the cylinder. A comparison of the axial displacement, 
velocity and acceleration will indicate the cylinder's accordion motion. Plots in the z- 
direction indicate the whipping motion of the cylinder. These time history plots will be 
presented later in this chapter to compare the discrete model's response to the simplified 
model's dynamic response. 
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lib) State time of 1.0 msec 




11c) State time of 1.5 msec 

Figure 11a, b,c. Deformation Results of Unstiffened ONR Cylinder Model 




lid) State time of 2.0 msec 




lie) State time of 2.5 msec 




Ilf) State time of 3.0 msec 



Figure lld,e,f. Deformation Results For Unstiffened ONR Cylinder Model (cont.) 
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B. SOR BEAM MODEL 



The SOR beam element formulation is based on approximating the shell's curvature 
with surfaces of revolution. This type of approximation has extensive applications for 
modeling containers, tanks, domes and has been done for a long period of time [Ref. 9]. 
A small displacement of a point on the surface of revolution can be resolved into two 
components. One of the components is in the direction tangent to the curvature of the 
surface about its SOR axis. The other component is in the direction of the normal to the 
beam surface. The displacement about the surfaces of revolution are estimated by applying 
a sine and cosine harmonic function for approximating the response on the surfaces of 
revolution. The harmonic function is then incorporated into the double asymptotic 
approximation procedure for estimating the fluid structure interaction in an infinite medium. 
As previously noted, DeRuntz [Ref. 6] developed an application to USA codes that assigns 
a harmonic function to each of the SOR elements or branches. The harmonic function will 
describe the SOR fluid degrees of freedom and more closely predict the responses to the 
underwater explosive loading. During this study, the beam was modeled with separate sine 
and cosine functions, however it was proven that a more accurate numerical result was 
obtained when both sine and cosine functions were used for modeling the beam with SOR 
elements. The net result was that the unstiffened and stiffened cylinder could be effectively 
modeled as a beam using the SOR element technique. 

The SOR beam model was modeled with the same number of finite elements (45) that 
was used to model the discrete cylinder along its axis. The model is shown in Figure 12 and 
was subjected to explosive loading. The beam had the same material properties as the 
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cylinder's shell and its cross-sectional area and moments of inertia were computed using the 
cylinder's dimensions. That is, the cross-section area of the SOR beam is 

A - \id] - df) 

4 

where d 0 and dj are the outside and inside diameters, respectively. Similarly the second 
moment of inertia and polar moment of inertia of the cross-section are 

1 = ^ d ° ~ ^ 

64 

and 

J - - df) 

x 

y 




Figure 12. SOR Beam Model With 45 Elements 
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The global deformation of the beam model was plotted for state times 0.5 msec to 3.0 
msec and is shown in Figure 13. The plots are magnified by a factor of 50 for enhanced 
visual recognition. There is an appearance of some crumpling or wrinkling at a few peaks 
with an increase in time. The deformation gradient at the ends are not as distinct as with the 
cylinder, however the overall shape and curvature is in excellent agreement with the discrete 
model. The displacement, velocity, and acceleration time histories were also generated and 
will be discussed later in the chapter. 
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13a) State time of 0.5 msec 





13b) State time of 1.0 msec 






13c) State time of 1.5 msec 






13d) State time of 2.0 msec 






13e) State time of 2.5 msec 






13Q State time of 3.0 msec 





Figure 13a,b,c,d,e,f. Deformation Plots for SOR Beam Model with 

50X magnification factor 
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c. 



INTEGRATED SOR BEAM/SHELL MODEL 



A third study was undertaken to simplify the original unstiffened cylinder by 
integrating SOR beam elements with shell elements. This type of model would not only 
reduce the size of the discrete model but would now be sufficiently small enough to allow 
for additional internal structures to be added to the shell section of the beam/shell model. 
The same harmonic functions for the SOR elements were used and the shell elements were 
generated in the same manner as in the discrete model. The success of this model solely 
depended on developing the proper interface between SOR beam and shell elements that 
would produce accurate dynamic response to the explosions. The first consideration for the 
interface between both types of elements was using internal beams to connect SOR elements 
to shell elements. This consideration proved to be totally ineffective because there was 
entirely too much local deformation at the beam/shell interface. 

The next type of interface explored was the use of a dry endplate with thickness equal 
to the thickness of the shell. More favorable dynamic responses were obtained, however 
there were sizeable differences in the axial motions of the model. This type of model is 
represented in Figure 14. The model is shown as a full model, even though a half model was 
used to reduce computational time and take advantage of symmetry conditions. 

1. Variation of Material Properties for Beam/Shell Interface 
A sensitivity study was undertaken by changing the material property of the dry 
endplate being used as the interface. Both the thickness of the plate and Young's modulus 
were adjusted in order to observe changes in the deformation, both axially and in the 
transverse or z-direction. The transverse or whipping motion of the model was not affected 
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X 





Figure 14. Integrated SOR Beam/Shell Model with Dry Endplate Used as 
Interface 

by large changes in the material properties for the endplate. The axially direction, however 
was affected by changes in either endplate thickness or Young's modulus. 

It was noted that the dynamic responses were more sensitive to changes in endplate 
thickness. Only when Young's modulus was adjusted in a large amount, there was a 
noticeable change in the dynamic response of the entire model. After performing numerous 
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numerical studies using the dry endplate, it was determined that the dry endplate produced 
favorable results for the modeling the integrated SOR beam/shell model. Furthermore, the 
endplate's material properties can be approximated by the following formula: 

3. P* a 4 * 

< o 

4* E* t 3 

where: 

P= max incident pressure(psi) 
a= radius of plate 
E= Young's modulus 
t= thickness of plate 

6 = tolerance of axial displacement for the endplate 

The formula was simplified from the exact solution for the deflection of a circular plate 
loaded at the center [Ref. 9]. Results from application of the formula have produced an 
integrated model that has very similar whipping motion responses and slight differences in 
axial direction responses. The reasons for slight differences in axial direction responses can 
be attributed to approximating the shell with SOR elements. Although SOR elements 
produce good whipping motions, the axial motions are slightly off because the shell elements 
not only bend axially but also bend in the radial direction causing differences in the overall 
result. The SOR beam model did not adjust locally for the mass of the endplates which also 
causes differences in the axial direction. Locally smearing the properties of the endplates 
will be examined in the discussion on stiffened beam/shell models. 

Figure 15 shows a comparison of the SOR/shell model with a 3/16" thick dry 
endplate. Figure 15a shows the model's global deformation at 0.5 msec after the shock has 
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hit the structure and before the interface endplate's material properties were adjusted. 
Figure 15b shows the global deformation after the endplate's material property was adjusted. 
Notice that the unadjusted model has a noticeable amount of axial deformation at the 
endplate as compared to the adjusted model which shows very little axial deformation at 
endplate. The tolerance for the adjusted model was set at .01 inch. Using a small tolerance 
ensures that there will be very little deflection of the dry endplate which improves the overall 
global deformation of the model. Also, the small tolerance ensures a smooth transition of 
deformation at the interface. There is some pinching of the shell caused by the rigidity of the 
dry endplate, however it will be seen later in this chapter that the overall global deformation 
and dynamic responses at the beam/shell interface are in good agreement with the discrete 
model. This type of model also provides for a smooth transition of the dynamic responses 
from beam to shell at the interface region. 



29 



15a) 3/16" Thick Endplate- No Adjustment to Material Properties 
50 X Magnification Factor , state time of 0.5 msec 




15b) 3/16" Thick Endplate, Material Property Adjusted, 
50 X Magnification Factor, state time 0.5 msec 



Figure 15a, b. Global Deformation Plots for SOR Beam/Shell Model: 
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The beam/shell model with adjusted material properties provides very similar 
deformation results when compared to the discrete model, however there is one disadvantage 
to this type of modeling. Adjusting the elastic modulus causes a longer computation time 
due to the increased stiffness in each element. The time step size used in the explicit finite 
element computer program is computed by the following formula: 



A t = 




where: 

S= safety factor (use number less than 1) 

L=length of smallest element in any direction 
E= Young's modulus 
p= material density 

By selecting a larger Young's modulus, the overall time step is decreased which increased 
the overall computational time. 

t 

After formulating the proper way to integrate SOR beam/shell elements, a study was 
done to find a dry endplate that would produce accurate results but reduce the computational 
time. A thicker endplate was first examined since it would help reduce the ratio of Young's 
modulus to material density. The Young's modulus for the 1/2" thick endplate was 
calculated using the endplate formula. A 1/2" thick dry endplate was selected for the 
integrated SOR/shell model since this produced an adjusted Young' modulus of about 100 
times less than the 3/16" thick endplate's Young's modulus. Both models were subjected to 
the same explosive testing with the charge location at the symmetric point on the shell. 
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Figures 16-18 compare the axial displacements, velocities and accelerations between the 
3/16" and 1/2" thick endplate models out to 3.0 msec. Since only elastic response was 
examined, the responses out to 3.0 msec produces the most effective time to compare the 
results. There is exceptional agreement with all three plots, however there is a slight 
difference in magnitude and frequency which can be accounted for by different masses of the 
endplates. The 1/2" thick plate produced slightly larger displacements in the axially 
direction. Figures 19-21 compares the z-direction or whipping motion responses of the two 
models at the shock point on both models. There is very tight agreement between both 
models for the displacement, velocity, and acceleration. This clearly suggests that the 
estimation of the material property of the endplates is a very effective tool for integrating the 
SOR Beam/Shell model. Even though a thicker endplate was used the overall time step 
was (9.027xl0‘ 8 ) which is 30 times less than the original time step of (2.5xl0 -6 ). Further 
optimization of this problem would produce a dry endplate of about 1.75" that would have 
a time step of about 80% of the original time step. Because DYNA3D code allows a material 
choice option for modeling a material as a rigid body, the time step problem can be 
eliminated. This material option will now be discussed. 
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displacement, (in.) 



Axial-Displacement For Two Beam/Shell Adjusted Models 




Figure 16. Axial Displacement of Two Integrated SOR Beam/Shell Models 
Both With Material Property Adjustment to Endplates 
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Axial- Velocity For Two Beam/Shell Adjusted Models 




Figure 17. Axial Velocity Response for Two Integrated SOR Beam/Shell Models 

X10 5 Axial-Acceleration For Two Beam/Shell Adjusted Models 




Figure 18. Axial Acceleration Response for Two Integrated SOR Beam/Shell Models 
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Z-Displacement For Two Beam/Shell Adjusted Models 




Figure 19. Whipping Displacement Response for Two Integrated Models 
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Z-Velocity For Two Beam Shell Adjusted Models 




Figure 20. Whipping Velocity Response for Two Integrated Models 




Figure 21. Whipping Acceleration Response for Two Integrated Models 
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The rigid body material option on DYNA3D is built into the computer code and 
only requires three material property inputs (E, p, v ). The thickness of the endplate was 
first selected at 3/16" and then varied to see how well this material option would predict the 
dynamic responses. It is recommended that the thickness of the endplate be very close to 
the thickness of the cylindrical shell or too much local deformation will occur at the interface 
area due to the increased mass distribution of the thicker endplate. The structure created 
with this material option will then behave as a rigid body dynamically when subjected to any 
type of loading. The SOR beam/shell model was redesigned using this material option and 
tested. Figure 22 shows the global deformation plots for the model using the rigid body 
material option. As with the previous modeling of the beam/shell interface, the rigid body 
model produces pinching at the endplate interface. Figure 23 shows a comparison plot 
between the adjusted 3/16" thick endplate model and the beam/shell model with the rigid 
body endplate at 0.5 msec. Both models were in excellent agreement with each other. The 
big advantage of using this material option is that the original values of Young's modulus and 
density can be used which in return does not affect the original time step. This significantly 
reduces the overall computational time of the model and makes it a viable choice for 
simplifying the cylindrical shell. 
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Figure 22b) State time of 1.0 msec 




Figure 22c) State time of 1.5 msec 

Figure 22a, b,c. Deformation Results for Integrated SOR Beam/Shell Model Using 
Rigid Body Material Option with DYNA3D 
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Figure 22d) State time of 2.0 msec 




Figure 22e) State time of 2.5 msec 




Figure 22f) State time of 3.0 msec 

Figure 22d,e,f. Deformation Results for Integrated SOR Beam/Shell Model (cont.) 
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Figure 23a) Deformation Plot of 3/16" Thick Endplate for SOR Beam/shell Model 




Figure 23b) Deformation Results of Rigid Body Endplate for SOR Beam/shell Model 



Figure 23 a,b. Deformation Results for 3/16" Thick Endplate SOR Beam/shell Model 
Versus Rigid Body Endplate SOR Beam/shell Model of Unstiffened 
ONR Cylinder 
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D. COMPARISON OF DISCRETE TO SIMPLIFIED MODELS 



This section will compare the results of the explosive loading on the discrete 
unstiffened model to those of the simplified SOR beam and integrated SOR beam/shell 
models. Comparisons will be done by plotting the displacement, velocity and acceleration 
in both the axial and transverse(z) directions for all three models at selected points. The z- 
direction or whipping motion response will be examined at the center point nearest to the 
charge. The axial or accordion motion responses will be monitored at the left endplate for 
all three models. 

For the cylinder, the center point nearest to the charge (so-called side shell node) and 
at the top of the cylinder at the center (called top shell node) and center point away from 
charge (called back side node) could be selected for comparison. Because the cylinder has 
both global and local deformations at the three nodal points, the locations are expected to 
have different responses. Depending on which motion is dominant, the difference may be 
large or small. If the global response is dominant over the local response, the three points 
will have a similar response. Otherwise, the three points will have quite different responses. 
In the z-direction, the global response dominanted over the local response. This can be 
clearly seen in Figure 24 which compares the three shell nodal points to the same beam nodal 
point on the simplified SOR beam model. Notice how the responses are very similar and that 
the top and back side shell nodes have a slight delay before the deformation occurs. This is 
due to the shock wave progressing through the cylinder where it hits the front side shell node 
first. Except for the delay, the responses are quite similar in shape and magnitude. This 
study will use the front side shell node only for comparison purposes. 
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Z-Displacement 




Figure 24. Whipping Displacement Plot Showing Comparison of Shell Node Response 
for Unstiffened Shell versus SOR Beam Model 
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The axial displacement, velocity, and acceleration time history plots for the discrete 
and two simplified models are provided in Figures 25-27. The beam/shell model was 
modeled with the rigid body material option using a 3/16" thick endplate for more accurate 
results and a quicker computational time. There is however noticeable differences between 
all three models for all three responses. An explanation for the differences is that the SOR 
beam model will not have the same local nodal point response as the same point on the 
cylinder. The local response in the axial direction is more dominant than the global response. 
The beam/shell model and the beam model did not account locally for the mass of the 
endplate which causes some error in the axial direction. If the magnitude of the local 
responses are examined, there is about a 45% difference between shell and beam model 
displacements. Figure 28 shows the global deformation plot at 0.5 msec for all three models. 
The overall global deformation is in good agreement except for local deformation at the 
beam/shell interface. 

Figures 29.-3 1 provide the whipping responses for all three models at the location of 
the closest point to the charge. These plots show excellent agreement between all three 
models with the error less than 10%. The SOR approximation predicts responses better in 
the transverse direction than in the axial direction. 
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0.05 



Axial-Displacement for All Three ONR Models 




Figure 25. Axial Displacement for All Three Models 
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Axial-Velocity for All Three ONR Models 




Figure 26. Axial Velocity for All Three Models 

vIO 6 Axial-Acceleration for All Three ONR Models 




Figure 27. Axial Acceleration Plot for All Three Models 
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28a) Global Deformation for Unstiffened Cylinder 




28b) Global Deformation for SOR Beam Model 




28c) Global Deformation Result For SOR Beam/Shell Model 



Figure 28a, b,c. Global Deformation Plots for All Three Models at 0.5msec 
with 50X Magnification Factor 
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Z-Dispiacement for All Three ONR Models 




Figure 29. Whipping Displacement for All Three Models 
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Z- Velocity for All Three ONR Models 




Figure 30. Whipping Velocity Plot for All Three Models 



Z-Acceleration for All Three ONR Models 




Figure 31. Whipping Accelearation Plot for All Three Models 
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And finally, Table 2 is provided to compare the size of the model with the overall 
computational time for each model. The SOR beam/shell model data is for the rigid body 
material option for modeling the beam/shell interface. Clearly, one must make a decision 
as to whether to use SOR beam elements or original shell elements. Clearly the accuracy lies 
in using shell elements, however if a section of the cylinder model requires isolation for 
detailed internal modeling, then the beam/shell model can be used to simplify the design and 
to reduce computational time. The next part of the study will examine modeling the ONR 
stiffnened cylinder. 





Unstiffened ONR 
Cylinder 


SOR Beam Model 


Integrated 
Beam/shell Model 


Number of 
Nodes 


1140 


47 


534 


Number of Wet 
Elements 


1080 


45 


360 


Computational 
Time (sec) 


1968 


6.44 


742 



Table 2. Comparison Data For all Three Models of Unstiffened ONR Cylinder 
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V. MODELING OF STIFFENED ONR CYLINDER 



In this section the simplification of the ONR cylinder is investigated. The original 
ONR cylinder with stiffeners will be modeled first. Next, the stiffener smearing technique 
will be used to model the cylinder with orthotropic shell elements. The SOR beam model 
will be generated next and finally, the integrated SOR beam/shell model will be generated. 
The original material used for all four models is as specified in Table 1. Only the smearing 
values will be different. This section will conclude with a comparison of the displacement, 
velocity and acceleration responses for the models subjected to both a symmetric and off- 
center shock wave. 

A. DISCRETE STIFFENED ONR CYLINDER 

The three-dimensional model is shown in Figure 32. As with the unstiffened model, 
a full model is depicted even though a half model was analyzed due to symmetry conditions. 
Both the unstiffened and stiffened models contain 1080 wet shell elements, however with 
the stiffeners, the total amount of elements for the stiffened model is 2064. The term wet 
means exposed to the shock wave pressure. This type of model involves longer 
computational time to complete the computation of the response to the shock wave due to 
the increased number of wet elements. It is because of both the computational time and 
design effort that a simpler model is needed. 

The stiffened or discretized model was subjected to both symmetric and off-center 
shock waves. As with the unstiffened model, there was no global deformation observed in 
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Figure 32. ONR Full Cylinder Model with Internal Stiffeners 

the x-axis (see Figure 32). The global deformation responses are provided in Figure 33 at 
discrete states of time. The plots are again magnified at a scale factor of 50 for viewer 
enhancement. The deformation analysis indicates that there is a large local deformation 
gradient at the locations of the deep frames and endplates. The gradient smoothed out as the 
shock wave progresses through the cylinder and by 3 msec, there is very little pinching effect 
caused by the deep frames. Displacement, velocity, and acceleration plots will be presented 
later in this section. 
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Figure 33a) State time of 0.5 msec 




Figure 33b) State time of 1.0 msec 




Figure 33c) State time of 1.5 msec 

Figure 33a, b,c. Deformation Results of Stiffened ONR Cylinder Model 
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Figure 33d) State time of 2.0 msec 




Figure 33e) State time of 2.5 msec 




Figure 33f) State time of 3.0 msec 

Figure 33d, e,f Deformation Results of Stiffened ONR Cylinder (cont.) 
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B. ORTHOTROPICALLY SMEARED ONR CYLINDER MODEL 



The first simplification model of the stiffened ONR cylinder investigated was the 
orthotropically smeared ONR cylinder. This type of simplification attempts to properly 
model the cylinder as one unit with homogenized properties. Once the effective or smeared 
properties are computed, the cylinder can be subjected to the underwater explosion and the 
results analyzed. The DYNA3D code allows the user to select an orthotropic elastic 
material. The material selection requires 9 properties in order to complete the material 
matrix. The six properties would represent the smeared properties of the original ONR 
cylinder. The first three properties were the Poisson's ratio for all three material axes. The 
material axis corresponds to the axis shown in Figure 1. All three of these values were 
assumed to be the same value of the original material or 0.3. 

Next the values for the Young's modulus were developed. The value for smeared 
Young's modulus in the axial direction (y-direction) is simply the discrete value for the 
original material since ring stiffeners do not contribute much in the longitudinal direction. 
To determine the smeared elastic modulus values in the radial (x or z-direction) and 
tangential direction (x or z-direction), the assumption was made that the radial stiffness of 
the shell is significantly less important than the tangential stiffness. The tangential stiffness 
in other words directly affects the bending stiffness of the cylinder. This assumption is also 
stated by Timoshenko and Woinowsky-Kreiger [Ref. 9]. Therefore the smeared value in 
the tangential direction would be set equal to the radial direction smeared value ( E e . = E, ). 
To calculate this smearing value for the radial and tangential direction, the 2-dimensional 
static smearing technique developed in [Ref. 9] was used. This smearing value would be 
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used in the three-dimensional representation of the stiffened ONR cylinder. The smearing 
value for both directions was computed by the following formula [Ref. 9] : 



E' = E ♦ [ 



12(1- v 2 ) 



a 



where: 

E'= E e . = E, 

E= discrete or original Young's modulus for cylinder 

v= Poisson's ratio 

h= thickness of cylinder 

E s = Young's modulus for stiffeners 

Is = moment of inertia of the stiffener 

The remaining three values are for the elastic shear modulii for all three planes. The values 
for G xy ' and G^' were assumed equal to each other due to the same reasoning for the elastic 
modulus. The value for these two planes was found by using the following formula [Ref. 9]: 




E 

2(l + v) 



where 

E= elastic modulus of original cylinder 
v= Poisson's ratio for original cylinder 

The value for the elastic shear modulus G n ' was found using the following expression 
[Ref.9]:: 
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where 



G= original elastic shear modulus of cylinder 
G s = shear modulus of stiffeners 
a= spacing between stiffeners 

Besides ensuring that the proper stiffness smeared values are inserted into the 
smeared model, the effective mass distribution must also be accounted for. The effective 
mass distribution for the smeared model was calculated by using the formula: 




where: 

p' = smeared mass density value 

M = total mass including cylinder and stiffeners to be smeared 
V= volume of cylindrical shell without stiffeners 

Due to the large thickness of the deep frame stiffener, the deep frames were discretely 
modeled in the smeared ONR cylinder model. The model is presented in Figure 34 with 
the axis as shown. Table 3 summarizes the inputs for the orthotropic model. This model 
was then subjected to the underwater explosive testing for both symmetric and off-center 
shot geometries. The overall global deformation plots for the symmetric shot geometry is 
provided in Figure 35 at a magnified scale factor of 50. There is slight pinching around deep 
frames at the earlier state times but the deformation gradient around each of the deep frames 
is less than the original discrete model. As will be seen in futures ection in this chapter, the 
overall dynamic reponses of this model are in close agreement with the discrete model. 
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Figure 34. Orthotropically Smeared Model of ONR Stiffened Cylinder 
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Table 3. Material Properties for Orthotropic Smeared ONR Cylinder 
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The use of the orthotropic smeared technique produces results similar to the study on 
the flat plate where local smearing was used around the clamped edges. In both cases, local 
smearing reduced the deformation gradient and produced effective global deformation 
results. One of the disadvantages of using the orthotropic smeared model is more lengthy 
computation times due to explicit time integration scheme. Although the model was 
reduced from 2064 nodes to 1124 nodes, the amount of wet elements remains the same. 
Since smearing also involves increasing the Young's modulus, the overall critical time step 
is decreased from the original isotriopic model's time step. These factors cause a very 
lengthy computational time. 
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Figure 35a) State time of 0.5 msec 




Figure 35b) State time of 1.0 msec 




Figure 35c) State time of 1.5 msec 

Figure 35a, b,c. Deformation Results of Orthotropically Smeared ONR Cylinder Model 
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Figure 35d) State time of 2.0 msec 




Figure 35e) State time of 2.5 msec 




Figure 35f) State time of 3.0 msec 

Figure 35d,e,f. Deformation Results of Orthotropically Smeared ONR Cylinder (cont.) 
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C. SOR BEAM MODEL 

The SOR beam model for the stiffened ONR cylinder is developed in the same 
procedure as the unstiffened SOR beam model. However, the effective stiffness and mass 
distribution of the stiffened model must be incorporated into the SOR beam model. This will 
ensure that the beam model will have the equivalent fluid structure interaction as compared 
with the stiffened ONR cylinder. The same harmonic function will be used in computing the 
fluid mass matrix for the SOR beam model. Both symmetric and off-center shot geometry 
will be used for analyzing the model. 

There are 45 SOR beam elements used to model the cylinder along its longitudinal 
direction. In order to obtain the equivalent stiffness, the cross-sectional areas and moments 
of inertia were adjusted. The cross-section area of the stiffened SOR beam model is: 

A - j-(d 2 - d 2 ) 



where d 0 and dj are the outside and inside diameters respectively. This value is the same for 
both stiffened and unstiffened models. However, the second moment of inertia and polar 
moment of inertia of the stiffened cylinder differed from the unstiffened cylinder because 
the stiffeners were considered. The values were computed as follows: 



[^M, 4 - df)\ 

4 



cylinder 



[-7 (< - d% eepframe ♦ [l(d 4 0 - d 4 )] shallowframe 
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J = 2. / 



The mass distribution value was the same as the orthotropically smeared mass value. 
Figure 36 shows the SOR beam model of the stiffened ONR cylinder. The model was 
subjected to shock wave and the deformation plots were made and shown in Figure 37. The 
deformation has a similar shape to the discrete model, however closer examination revealed 
that the deformation was in fact less than the discrete model. This was due to not taking into 
account the local effects of the deep frames and endplates. Since the mass distribution and 
cross-sectional properties of the model only accounted for the shallow frames , the model 
doe not properly show the local mass density and equivalent stiffness of the endplates and 
deep frames. 
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Figure 36. SOR Beam Model of Stiffened ONR Cylinder 
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Figure 37a) State time of 0.5 msec 



Q 



Figure 37b) State time of 1.0 msec 



Figure 37c) State time of 1.5 msec 



Figure 37a, b,c. Deformation Results of SOR Beam Model of ONR Cylinder Model 
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Figure 37d) State time of 2.0 msec 
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Figure 37e) State time of 2.5 msec 
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Figure 37f) State time of 3.0 msec 

Figure 37 d,e,f. Deformation Results for SOR Beam Model of ONR Cylinder (cont.) 
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To correct the mass distribution and equivalent stiffness of the SOR beam model, 
the mass density of the deep frames and endplates were inputted at their respective element 
location in the DYNA3D program. The cross sectional properties were also adjusted 
locally for the deep frames and endplates. Additionally, the element length for the deep 
frame section was reduced to two inches which corresponds with the original cylinder 
dimensions. This allowed for a more accurate local smearing of both mass density and cross- 
sectional properties for the deep frames. The adjusted SOR Beam model now contains 49 
elements and is shown in Figure 38. The global deformation plots are presented in Figure 
39. The deformation plots show an improvement from the globally mass smeared 45 SOR 
beam element model. Figures 40-42 show the axial dynamic responses at the left endplate 
which support using local smearing of the mass density for deep frames and endplates. 
Detailed comparisons between discrete and SOR beam model will be presented later in this 
report. The computational time using the SOR beam model is approximately 1/20 of the 
discrete model’s computational time. 

x 
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Figure 38. SOR Beam Model With Local Mass Density Smearing 
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Figure 39a) State time of 0.5 msec 




Figure 39b) State time of 1.0 msec 
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Figure 39c) State time of 1.56 msec 



Figure 39a, b,c. Deformation Results For SOR Beam Model With Local Mass Density 
Smearing of Endplates and Deep Frames 
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Figure 39d) State time of 2.0 msec 




Figure 39e) State time of 2.5 msec 




Figure 39f) State time of 3.0 msec 



Figure 39d,e,f. Deformation Results for SOR Beam Model With Local Mass Density 
Smearing of Deep Frames and Endplates (cont.) 
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Axial-Displacement for All Three ONR Models 




Figure 40. Axial Displacement Response Comparing Local to Global Mass Density 
Smearing of SOR Beam Model 
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Figure 41. 



Figure 42. 



Axial-Velocity for All Three ONR Models 




Axial Velocity Response Comparing Local to Global Mass Density 
Smearing of SOR Beam Model 
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Axial-Acceleration for Al! Three ONR Models 




SOR Beam-Local Smear 
* ONR Orig Model 
SOR Beam-Global. Smear 
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Time, ( milli seconds) 



Axial Acceleration Response Comparing Local to Global Mass Density 
Smearing of SOR Beam Model 
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D. INTEGRATED SOR BEAM/SHELL MODEL 



The integrated SOR beam/shell model incorporates all simplification techniques 
discussed in this paper. The interface will be modeled using the rigid body material option 
and local smearing will be done for the beam elements to take into account the mass 
distribution of the deep frames and endplates. Since it was desired to examine the dynamic 
responses at one of the deep frames, a beam/shell model was produced with fewer shell 
elements as shown in Figure 43. The model contains a total of 450 shell elements and 29 
SOR beam elements. The shell section has the shallow frames and deep frames discretely 
modeled to produce the fastest computational time. The beam section has the cross-section 
and second moment of inertia values adjusted to account for the shallow and deep frame 
stiffeners. This model will be tested with a charge located off-center as shown in Figure 
44. The symmetric test on beam/shell model would have only produced a larger model and 
hence it was not tested. 

x 
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Figure 44. Offcenter Shot Geometry for Integrated Beam/shell Model 

The global deformation plots are provided in Figure 45. The deformation plot shows 
local pinching at deep frames and endplates which is consistent with the discrete model. The 
rigid body model of the dry endplate produces larger local deformation at the endplate. The 
overall global deformation appears to be in good agreement with the other three models. 

The computational time for the model was about 1/3 of the discrete model which 
clearly supports using this type of model. It is highly recommended to use as many SOR 
Beam elements as possible in the model design to reduce the computational time. Detailed 
analysis at various locations on the beam/shell model will be provided later in this chapter. 
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Figure 45a) State time of 0.5 msec 




Figure 45b) State time of 1.0 msec 




Figure 45c) State time of 1.5 msec 



Figure 45a, b,c. Deformation Results of Integrated SOR Beam/Shell Model of ONR 
Cylinder Model 
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Figure 45d) State time of 2.0 msec 




Figure 45e) State time of 2.5 msec 




Figure 45f) State time of 3.0 msec 

Figure 45d,e,f Deformation Results for Integrated SOR Beam/Shell Model (cont.) 
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E. COMPARISON OF DISCRETE TO SIMPLIFIED MODELS 



This section will compare the dynamic responses between the fully discretized 
stiffened ONR cylinder to the simplified models. The analysis will compare both symmetric 
and off-center geometry for the location of the charge. The displacement, velocity, and 
acceleration responses will be displayed and discussed. 

1. Comparison Using Symmetric Explosive Testing 
The geometry for symmetric testing is the same as Figure 2. The aim of this testing 
was to simply compare the responses between the orthotropically smeared, beam, and 
discrete models of the stiffened ONR cylinder. The beam/shell model was not included since 
more detailed data was taken using off-center shot geometry. The results will aim to 
compare the responses at a deep frame location, endplates and at the shot location. 

Figure 46 compares the overall global longitudinal deformation plots for all four 
models at state time 0.5 msec and at magnified scale factor of 50. All plots are taken at the 
same coordinate axis position. There is good agreement except for areas of local 
deformation at the endplates and deep frames. The deformation plots show that both the 
smeared stiffener technique and SOR beam element modeling are effective in producing 
similar global responses. However, it is recommended that use of large numbers of 
orthotropic shell elements should be avoided since longer computational times occur when 
using these types of elements due to the decreased critical time step. 
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Figure 46a) ONR Stiffened Cylinder at 0.5 msec 




Figure 46b) Orthotropically Smeared ONR Cylinder at 0.5 msec 



Figure 46c) SOR Beam Model of ONR Cylinder at 0.5 msec 



Figure 46a, b,c Deformation Results for Discrete vs. Simplified Models for ONR 
Stiffened Cylinder with Symmetric Underwater Explosion Testing 



Next, the axial displacement, velocity, and acceleration plots were made at the left 
and right endplate areas to check global motion at the ends. Figures 47-49 shows the 
displacement, acceleration of all four models at the left endplate. There is a shift in the 
frequency between the SOR beam and discrete model, however the beam's peak magnitude 
is very close to the shell's first oscillation peak magnitude. There is significantly less 
vibration and acceleration experienced in the beam model as compared to the shell model 
which is expected. The orthotropic model compares favorably to the discrete model in shape 
but is slightly off in magnitude. 



Axial-Displacement for All Three ONR Models 




Figure 47, Axial Displacement Comparison at Left Endplate 
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Axial- Velocity for All Three ONR Models 




Figure 48. Axial Velocity Comparison At Left Endplate 




Figure 49. Axial Acceleration Comparison at Left Endplate 
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Figure 50-52 shows the whipping responses for all three models at the shock point 
(nearest structure point from the charge). The orthotropic shell model and discrete model 
are very similar in magnitude and shape and have an average error of less than 15%. There 
is less velocity and acceleration in the SOR Beam model because the local smearing 
represents the average stiffness which’ produce fewer oscillations than the original model. 
The cylindrical models have radial deformation which produce more velocity and 
acceleration responses as the shock wave progresses through the shell. 



Z-Direction Displacement for all Three ONR Models 




Figure 50. Whipping Displacement Comparison at Shock Point 
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Figure 51. Whipping Velocity Comparison at Shock Point of All Three Models 



Z-Direction Acceleration for All Three ONR Models 




Figure 52. Whipping Acceleration Comparison at Shock Point For All Three Models 
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And finally, the dynamic responses are analyzed at a deep frame on the three models. This 
location is 75 inches from the center of the cylinder (see Figure 1). Only the whipping 
motion responses will be presented, however there was good agreement in the axial direction 
at the deep frame. There was reduced velocity and acceleration in axial direction which was 
consistent with the endplate responses. Figures 53-55 show the whipping responses at the 
deep frame. The results show excellent agreement between the discrete and SOR beam 
model with only slight differences in velocity and acceleartion for reasons previously noted. 



Z-Direction Displacement for all Three ONR Models 




Figure 53. Whipping Displacement Comparison at Deep Frame 
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Z-Direction Velocity for All Three ONR Models 




Figure 54. Whipping Velocity Comparison at Deep Frame Location 




Figure 55. Whipping Acceleration Comparison At Deep Frame 
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Overall, the dynamic responses can support using either SOR beam elements, or 
orthotropic shell elements. There are tradeoffs to be considered for choosing which 
technique to use. The orthotropic shell element produces more accurate responses than the 
SOR beam element, however the SOR beam element's computational time is about 20 times 
faster than an orthotropic shell elements. When only 49 elements are used for SOR beam 
model vice 1080 for the orthotropic shell model, there is a big reduction in computational 
time, but with the reduced accuracy. So both techniques should be considered for 
analyzing, then if time is the more important factor in the solution, use of SOR beam 
elements are recommended. Table 4 provides a comparison based on size of model and 
overall computational time for all three models subjected to symmetrci shock testing. 





Stiffened ONR 
Cylinder 


SOR Beam Model 


Orthotropic 
Smeared Model 


Number of 
Nodes 


2064 


51 


1224 


Number of Wet 
Elements 


1080 


49 


1080 


H Computational 
f Time (sec) | 


2070 


95 


12200 



Table 4 . Comparison Data For all Three Models of Stiffened ONR Cylinder 
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2. Comparison Using Off-center Explosive Testing 

This section will compare the dynamic responses between the fully discretized 
stiffened ONR cylinder to the SOR Beam and integrated beam/shell model for an off-center 
shot geometry. The ONR stiffened model was subjected to the underwater explosion at a 
location 25 inches to the right of the center of the cylinder (See Figures 1 &44). The SOR 
beam model was the same model used in the symmetric analysis but subjected to the off- 
center shock wave. Finally, the integrated beam/shell model was subjected to the same 
testing. The integrated beam/shell model used the rigid body material option for the endplate 
and locally smeared the mass distribution and cross-sectional properties of the deep frames 
and endplates. The orthotropic shell model was not considered due to its lengthy 
computational time. 

Figure 56 is a global deformation plot of all three models at 0.5 msec. The models 
appear to be in good agreement except for some local pinching at the dry endplate on the 
integrated beam/shell model. There is however some difference in shape between the SOR 
beam model and discrete model which is mainly due to applying the average stiffness and 
mass distribution at the deep frames and endplates. This approximation produces similar 
results but will not exactly match the discrete model. 

Since the main objective in this study was to simplify the original model, it is 
important to select an area on the shell section of the integrated beam/shell model and 
compare the responses between discrete and the simplified model. If the responses are in 
good agreement at this selected area, then further internal structures could be placed inside 
the shell section of the beam/shell model. Thus, only a small portion of the original model 
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Figure 56a) ONR Stiffened Cylinder at 0.5 msec 






Figure 56b) SOR Beam Model of ONR Cylinder at 0.5 msec 




Figure 56c) 



Figure 56a, b,c. 



Integrated SOR Beam/Shell Model of ONR Cylinder at 0.5 msec 

Global Deformation Plots - Off-center Explosion Testing for Stiffened 
ONR Cylinder 
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would need to be modeled discretely, and the rest of the model could be modeled with SOR 
beam elements. The selected area for comparison purposes was the deep frames which is 
located 75 inches right of the center of the original cylinder and is discretely modeled inside 
the shell section of the beam/shell model. Both axial and whipping responses were 
analyzed. Figures 57-59 show the axial responses and Figure 60-62 show the whipping 
responses at the deep frame. Since the integrated beam/shell model used SOR beam 
elements for 60% of the longitudinal length of the original cylinder, slightly higher 
inaccuracies were experienced. The integrated beam/shell model responses are in closer 
agreement to the discrete model. The whipping responses are very similar in shape and 
only differ in magnitude due to the use of SOR beam elements. Table 5 is provided for 
comparing the size and computational time for all three models. 





Stiffened ONR 
Cylinder 


SOR Beam Model 


Integrated 
Beam/shell Model 


Number of 
Nodes 


2064 


51 


942 


Number of Wet 
Elements 


1080 


49 


479 


Computational 
Time (sec) 


2070 


95 


719 



Table 5. Comparison Data For all Three Models of Stiffened ONR Cylinder 
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displacement, (in.) 



0.14 



Axial-Displacement for All Three ONR Models 




Figure 57. Axial Displacement Comparison At Deep Frame Location 
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Axial-Velocity for All Three ONR Models 




Figure 58. Axial Velocity Comparison At Deep Frame Location 

v ^ n 6 Axial-Acceleration for All Three ONR Models 




Figure 59. Axial Acceleration Comparison At Deep Frame Location 
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displacement, (in.) 



2-Direction Displacement for all Three ONR Models 




Figure 60. Whipping Motion Comparison At Deep Frame Location 



Z-Direction Velocity for All Three ONR Models 




Figure 61. Whipping Motion Comparison At Deep Frame Location 



Z-Direction Acceleration for All Three ONR Models 




Figure 62. Whipping Motion Comparison At Deep Frame Location 
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Overall, the responses favor using an integrated beam/shell model if additional 
internal structures require modeling. This model provides both quicker computational time 
and required accuracy of dynamic responses. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



Simplified finite element modeling of a stiffened structure is an effective means by 
which both the design effort and the computational time can be reduced. When the discrete 
and simplified models are subjected to underwater explosions and closely examined, the 
results show that either smearing orthotropically or using the SOR beam element produces 
similar global results. When examining the two models at very discrete points, the results 
are similar but differ. The orthotropically smeared cylindrical model produces more accurate 
comparisons, but has a very lengthy computational time. There is more inaccuracy using 
SOR beam element modeling, however the computational time and design effort is 
significantly reduced. 

Incorporating the advantages of both techniques into an integrated beam/shell model 
provides the best means of simplifying the discrete model. The results are more favorable 
and the computational time is still significantly reduced. The interface between the beam 
and shell section can be modeled as either a dry endplate with smeared properties or by using 
the rigid body material option within DYNA3D code. It is recommended that the rigid body 
material option be used since the critical time step is not changed when this material is 
selected for modeling the interface. The use of the beam/shell model allows further internal 
structures to be modeled without dramatically increasing both design effort and 
computational time for numerical solution. If a quick but less accurate solution is needed 
for a particular problem, then the SOR beam model is recommended. 
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Since no physical experimental data was used in this study, it is highly recommended 
that a scaled down experimental model be built and tested against the numerical data 
collected in this study. The present numerical analysis involved using an explicit time 
integration scheme adopted in DYNA3D code. An additional comparison of the results 
could be accomplished by using an unconditionally stable implicit time integration 
technique with another finite element program. 
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APPENDIX A 



INPUT FILES FOR SOR BEAM MODEL 
The following input files consist of ingrid which creates the geometric model 
for DYNA3D, flumas, augmat, and timint. The flumas, augmat and timint files 
provide input for the USA program. These files are for SOR Beam Model Only. 
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T T T T 


$ 


PRTGMT 


PRTTRN 


PRTAMF 


CALCAM 


T F F F 


$ 


EIGMAF 


TWODIM 


HAFMOD 


QUAMOD 


F F F T 


$ 


PCHCDS 


NASTAM 


STOMAS 


STOINV 


F F F F 


$ 


FRWTFL 


FRWTGE 


FRWTGR 


FRESUR 


F T F F 


$ 


RENUMB 


STOGMT 


ROTGEO 


ROTQUA 


F F F F 


$ 


PRTCOE 


STRMAS 


SPHERE 


ROTSYM 


F F F F 


$ 


OCTMOD 


CAVFLU 


FRWTFV 


INTCAV 


DYNA 


$ 


MAINKY 








0 51 0 0 


$ 


NSTRC 


NSTRF 


NGEN 


NGENF 


10 0 


$ 


NBRA 


NCYL 


NCAV 




0 112 0 


$ 


NHAS 


NHAF 


NHAI 


NFUN 


51 


$ 


NSEG 








0 


$ 


NCIR 








9.356E-05 57600.0 


$ 


RHO 


CEE 






5 


$ 


NVEC 








1 50 16.8 16.8 3 


$ 


N1 


N2 


R1 R2 


NSET 



110. 16.8 1 
1 50 16.8 16.8 49 

50 50 16.8 0. 1 
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AUGMAT DATA FOR MODEL . BEAMONR 



strnam.dac 

P P P P 


beam . daa beam. geo beam.pre 


$ STRNAM FLUNAM GEONAM PRENAM 


r r r r 

C *T r* *T 




$ FRWTGE FRWTST FRWTFL PLNWAV 


r 1 r I 

D n n ^ 


$ 


FLUSKY 


DAAFRM 


SYMCON 


DOFTAB 


r t t I 

p p n n 


$ 


PRTGMT 


PRTTRN 


PRTSTF 


PRTAUG 


r r r r 
F 


$ 


MODTRN 


STRLCL 


INTWAT 


CFADYN 


DYNA 


$ 


MAINKY 








0 . 0 

^ 


$ 


NTYPDA 









r=d3dump 

TIME INTEGRATION DATA FOR MODEL BEAMONR 



beam.pre beam. 


pos 


$ 


PRENAM 


POSNAM 






beam. rst 




$ 


RESNAM 








F F F 




$ 


BUB PUL 


REFSEC 


FLUMEN 




T F F F 




$ 


EXP WAV 


SPLINE 


VARLIN 


PACKET 


F T F F 




$ 


HYPERB 


EXP LOS 


DOUBDC 


VELIMP 


1 




$ 


NCHARG 








0 . 0 




$ 


HYDPRE 








0.0 -25.3125 326.4 


$ 


XC YC 


ZC 






0.0 -25.3125 


16.813 


$ 


SX SY 


SZ 
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$ 


JPHIST 








1 . 0 




$ 


PNORM 








14 -0E-06 




$ 


DTHIST 








1 




$ 


CHGTYP 








60.0 25.8 10 


. 0 


$ 


WEIGHT 


SLANT 


CHGDEP 




1 0 




$ 


NT I NT 


NCHGAL 






0.0 2.00E-06 




$ 


STRTIM 


DELTIM 






5000 5000 




$ 


NS AVER 


NRESET 






0 0 0 0 




$ 


LOCBEG 


LOGRES 


LOCWRT 


NSTART 


F F 




$ 


FORWRT 


STBDA2 






F 




$ 


DISPLA 
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APPENDIX B 



INPUT FILES FOR SOR BEAM/SHELL MODEL 
The following input files are for the integrated SOR Beam/shell model. 



bsmodel .half 



9 


942 


0 


32 


860 


0 


0 


oo large 


0 0 


0 0 


0 














0 


0 e20 . 9 




450 33 






2 0 


450 0 


0 0 


0 0 


0 0 


0 0 


0 0 


0 0 


3 . 000E-03 


1 . 000E+0 3 


1 . 000E-05 




. OOOE+OO 


. OOOE + OO 


0 


. OOOE+OO 


0 0 


0 0 


0 2 


0 0 


1 0 


0 


0 0 


0 


0 0 


-OOOE+OO 














0 0 


0 














1 37 


. 3560E-04 


0 0 


.00 0 0E + 0 0 


0 . OOOOE 


+00 .OOOOE+OO 2 


0 0 


material type #3 (< 


elastic-plastic) 










2 . 900E+07 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE + OO 


.OOOE+OO 


-OOOE+OO 


.OOOE+OO 


3 .000E-01 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE + OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


2. 000E+10 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE + OO 


. OOOE + OO 


.OOOE+OO 


. OOOE+OO 


S . 020E+04 


. 000E+00 


.OOOE+OO 


-OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. 000E+00 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


- OOOE+OO 


.OOOE+OO 


. OOOE + OO 


.OOOE+OO 


. 000E+00 


. 000E+00 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE + OO 


.OOOE+OO 


section properties 

c 














1.880E-01 


1 . 880E-01 


1 . 880E-01 


1 . 880E-0 1 


. OOOE+OO 








2 17 


. 3560E-04 


0 0 


. 0000E+00 


0 .OOOOE 


+00 .OOOOE 


+ 

o 

o 

H* 


0 2 


material type # 1 (elastic) 












2 . 900E+07 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


3 . 000E-01 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. 000E+00 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE + OO 


. OOOE + OO 


. OOOE+OO 


.000E+00 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. 000E + 0 0 


.OOOE+OO 


• .OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE + OO 


.OOOE+OO 


. 000E+00 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


section properties 

n 














1 . 970E+0 1 


5 . 402E+03 


5 . 402E+03 


1.081E+04 


1 . 970E+01 


. OOOE + OO 






3 209 


.571OE-04 


0 0 


. OOOOE + O 0 


0 -OOOOE 


+00 .OOOOE 


+ 00 2 


0 C 


material type # 20 


(rigid body) 










2. 900E+07 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


3.000E-01 


•OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. 000E+00 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


-OOOE+OO 


. OOOE+OO 


. 000E+00 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. 000E+00 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. 000E+00 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE + OO 


.OOOE+OO 


. OOOE+OO 


section properties 

c 














1.880E-01 


1 .880E-01 


1 . 880E-01 


1 .880E-01 


. OOOE+OO 








4 37 


. 3560E-04 


0 0 


. OOOOE + OO 


0 .OOOOE 


+00 .OOOOE+OO 2 


0 0 


material type # 3 (elastic-plastic) 










2 . 900E + 0 7 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


3 . OOOE-Ol 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


-OOOE+OO 


.OOOE+OO 


.OOOE+OO 


2. OOOE+OS 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


-OOOE+OO 


. OOOE+OO 


S. 100E+04 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.000E+00 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


section properties 

c 














2 . S00E + 00 


2 . 500E+00 


2 . 500E+00 


2 . 500E+00 


. OOOE+OO 








5 37 


.3560E-04 


0 0 


. OOOOE+OO 


0 .OOOOE 


♦00 .OOOOE+OO 2 


0 o 


material type # 3 (elastic-plastic) 










2 . 900E+07 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE + OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


3 .000E-01 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


2 . 000E+05 


. OOOE+OO 


.OOOE+OO 


-OOOE+OO 


.OOOE+OO 


. OOOE + OO 


. OOOE+OO 


. OOOE+OO 


S . 100E + 04 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. 000E+00 


. OOOE+OO 


.OOOE+OO 


-OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 
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. 000E+00 


.000E+00 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


section properties 

c 












3 . 250E+00 


3 . 2S0E+00 


3 . 250E+00 


3 . 250E+00 


.OOOE+OO 








6 37 


. 3560E- 04 


0 0 


. OOOOE+OO 


0 .OOOOE+OO .OOOOE+OO 2 


0 0 


material type # 3 (elastic-plastic) 










2. 900E+07 


. 00 0E+ 0 0 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


3 . OOOE-Ol 


. 000E + 0 0 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


2 . 000E+05 


. 000E+00 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


5 . 100E+04 


.000E+00 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.000E+00 


.000E+00 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. 000E+00 


. 000E+00 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


section properties 














2 . 000E+00 


2 . 000E+00 


2 .OOOE+OO 


2. OOOE+OO 


. OOOE+OO 








7 37 


. 3560E-04 


0 0 


. OOOOE+OO 


0 .OOOOE+OO .OOOOE+OO 2 


0 C 


material type # 3 (elastic-plastic) 










2 . 900E+07 


. 000E+00 


. OOOE+OO 


. OOOE+OO 


. OOOE + OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OC 


3. 000E-01 


. 000E+00 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


2 . 000E+05 


. 000E + 0 0 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE + OO 


. OOOE+OO 


5 . 100E+04 


. 000E+00 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE + OO 


. OOOE+OO 


.000E+00 


. 000E+00 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.000E+00 


. 000E+00 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


section properties 














2 . 500E-01 


2 . 500E-01 


2 . 500E-01 


2.500E-01 


.OOOE+OO 








8 18 


. 2692E-04 


0 0 


. OOOOE+OO 


0 .OOOOE+OO .OOOOE+OO 1 


0 2 


material type # 1 (elastic) 












2 . 900E+07 


. 000E+00 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


. OOOE + OO 


. OOOE+OC 


3 . 000E-01 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OC 


. 000E+00 


. 000E+00 


. OOOE+OO 


. OOOE+OO 


. OOOE + OO 


. OOOE+OO 


. OOOE + OO 


. OOOE+OO 


. 000E+00 


. 000E+00 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. 000E+00 


. 000E+00 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.000E+00 


. 000E+00 


. OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


section properties 

n 














3 . 970E+02 


3 . 614E+03 


3 . 614E+03 


7. 227E+03 


3 . 970E+02 


. OOOE+OO 






9 11 


. 0470E-03 


0 0 


. OOOOE+OO 


0 .OOOOE+OO .OOOOE+OO 1 


0 2 


material type # 1 (elastic) 












2 . 9 00E+07 


. 000E+00 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


3 . 000E-01 


.OOOE+OO 


. OOOE+OO 


*. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE + OO 


. OOOE+OO 


. 000E+00 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. 000E+00 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. 000E+00 


. OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. 000E+00 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


. OOOE+OO 


.OOOE+OO 


. OOOE+OO 


.OOOE+OO 


section properties 

n 














8 . 880E+02 


6 . 275E+04 


6 . 275E+04 


1. 255E+05 


8 . 880E+02 


.OOOE+OO 







1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 



1 . 
0 . 
0. 
0 . 
0. 
0 . 
1 . 
0 . 
0 . 
0 . 
0 . 
0 . 



. 0000000000000E+00 -. 12656250000 00E+ 03 
.2630132913S89E+01 - . 12656250000 00E+03 
. 5195503234863 E+ 01 - . 12656250000 00E+03 
. 763 2 9431 5338 1E+ 01 - . 1265625000000E+03 
. 9882434 844 97 1E+ 01 - . 126S62SOOOOOOE+03 
. 11888 58 79S166E+02 - . 12656250 000 00E+03 
.000 00 00 000 000E+00 - .12093750000 00E+03 



.2630133152 00 8 E+ 01 
. 5195 504 188538 E+ 01 



-.1209375000000E+03 



-. 1681300163269E+02 
- .16606 0066223 1E+0 2 

- . 1599011611938E+02 

- .14 98049640 656E+02 

- . 1360200595856E+02 
-.1188858795166E+02 

- .16813 00163269E+02 

- . 16606008 52 96 6E + 0 2 



- .120937S0000 00E+03 - . 1599011707306E+02 



. 7632942199707E+01 - . 1209375000000E + 03 - . 1498049640656E+02 
. 9882434844971E+01 - . 1209375000000E + 03 - . 1360200691223E+02 
. 1188858795166 E+ 02 - . 1209375000000E+03 - . 118885879S166E+02 



5. 
0 . 
0. 
0. 
0 . 
0 . 
5. 
0 . 
0 . 
0 . 
0 . 
0 . 
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913 

914 

915 

916 

917 

918 

919 

920 

921 

922 

923 

924 

925 

926 

927 

928 

929 

930 

931 

932 

933 

934 

935 

936 

937 

938 

939 

940 

941 

942 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 



0. . 1259073257446E+02 
0. . 1386673641205E+02 
0. 14 8 012 952 804 6E+0 2 
0. . 1537139701843E+02 
0. . 1556300067902E+02 
0. . 153 713 9511 108E+02 
0. . 1480129432678E+02 
0. . 138667364 12 05E+02 

0. . 12S9073352814E+02 

1. .0000000000000E+00 
0. . 2434589624405E+01 
0. .4809231758118 E+ 01 
0. . 70 654 54 9 59 869E+01 
0. . 9147703170776E+01 

0 . . 11004703 52 173 E+ 02 

1. . 0000000000000E+00 
0. .2434589624405 E+ 01 
0. . 4809231758118E+01 
0. .7065454959869 E+ 01 
0. . 9147703170776E+01 
0. . 1100470352173E+02 
0. . 1259073257446E+02 
0. . 1386673641205E+02 
0. . 14 8012952 8046E+02 
0. . 1537139701843E+02 
0. .1556300067902 E-t-02 
0. . 1537139511108E+02 
0. . 14 80129432 678E+02 
0. . 138667364120SE+02 
0. . 1259073352814E+02 

ELEMENT CARDS 



- .1153125000000E + 03 

- .1153 1250000 0 0E+0 3 
- .1153125000000E+03 
- . 115312 5000000E+0 3 
-.1153125000000E+03 

- . 1153125000000E+03 
-.1153125000000E+03 
* -1153125000000E+03 

- . 1153 12 50 000 00E+ 03 
- . 1209375000000E+03 

1209375000000 E+03 

- . 120 9 37500000 0E+0 3 

- .1209375000000E+03 

- . 120 9 3 75 00 0 0 00 E+0 3 

- . 1209 3750 00 000E + 03 

- . 120 9 3750000 00 E+0 3 

- .1209375000000E+03 

- .1209 375000 000E+03 

- . 1209375000 00 0 E+03 

- . 120937 50 00000E+03 

- . 120937 50 00000E+ 0 3 

- . 1209 375000 000E + 03 

- . 1209375000 00 0E+03 

- .1209375000000E + 03 

- . 1209375000000E + 03 

- . 1209375000000E+0 3 

- . 120937 50 00000E + 03 

- . 1209375000 000E+03 
-.1209375000000E+03 
-.1209375000000E+03 

OR BEAM ELEMENTS -- 



9 


553 


554 


552 


2 


554 


555 


552 


2 


555 


556 


552 


2 


556 


557 


552 


2 


557 


558 


552 


2 


558 


559 


552 


2 


559 


560 


552 


2 


560 


561 


552 


2 


561 


562 


552 


2 


562 


563 


552 


8 


563 


564 


552 


2 


564 


565 


552 


2 


565 


566 


552 


2 


566 


567 


552 


2 


567 


568 


552 


2 


568 


569 


552 


2 


569 


570 


552 


2 


570 


571 


552 


2 


571 


572 


552 


2 


572 


573 


552 


2 


573 


574 


552 


8 


574 


575 


552 


2 


575 


576 


552 


2 


576 


577 


552 


2 


577 


578 


552 


2 


578 


579 


552 


2 


579 


580 


552 



. 9147703170776E + 01 
. 7065455436707E+01 
.4809232711792E+01 
.2434590816498E+01 
.OOOOOOOOOOOOOE+OO 
. 2434589624405E+01 
.480923 17 58118E+01 
-7065454959869E+01 
. 9147703170776E + 01 
.1556300067902E+02 
. 1537139511108E+02 
. 1480129432678E + 02 
. 138667364 120 5 E+ 0 2 
. 12590733 52 814 E+0 2 
. 11004 70352 173 E+ 0 2 
. 15563 000 679 02 E+0 2 
. 1537139511108E+02 
. 14801294 32678E+0 2 
. 1386673641205E+02 
. 1259073352814E+02 
. 1100470352173E+02 
.9147703170776E+01 
.7065455436707E+01 
.4809232711792E+01 
.2434590816498E+01 
.0000000000000E+00 
.2434589624405E+01 
.480923 17 58 118E+01 
.7065454959869E+01 
. 9147703170776E+01 



103 



2.500E-01 2 . 500E-01 2 

839 7 920 

2.500E-01 2 . 500E-01 2 

840 7 921 

2.500E-01 2 . 500E-01 2 

841 7 7 

2.500E-01 2 . 500E-01 2 

842 7 8 

2.500E-01 2 . 500E-01 2 

843 7 9 

2.500E-01 2.500E-01 2 

844 7 10 

2.500E-01 2 . 500E-01 2 

845 7 11 

2.500E-01 2 . 500E-01 2 

846 7 928 

2 . 500E-01 2 . 500E-01 2 

847 7 929 

2.500E-01 2 . 500E-01 2 

848 7 930 

2.500E-01 2.500E-01 2 

849 7 931 

2.500E-01 2.500E-01 2 

850 7 932 

2.500E-01 2 . 500E-01 2 

851 7 927 

2.500E-01 2 . 500E-01 2 

852 7 934 

2 . 500E- 0 1 2 . 500E-01 2 

853 7 935 

2.500E-01 2 . 500E-01 2 

854 7 936 

2.500E-01 2 . 500E- 0 1 2 

855 7 937 

2.500E-01 2.500E-01 2 

856 7 938 

2.S00E-01 2 . 500E-0 1 2 

857 7 939 

2.500E-01 2.500E-01 2 

858 7 940 

2.500E-01 2 . 500E-01 2 

859 7 941 

2.500E-01 2 . 500E-01 2 

860 7 942 

2.500E-01 2 . 500E-01 2 



500E-01 


2 . 500E-01 




921 


383 


364 


500E-01 


2 . 500E-01 




912 


132 


383 


500E-01 


2 . 500E-01 




922 


923 


8 


500E-01 


2 . 500E-01 




923 


924 


9 


50 OE - 0 1 


2 . 500E-01 




924 


925 


10 


500E-01 


2 . 500E-01 




925 


926 


11 


500E-01 


2 . 500E-01 




926 


927 


12 


500E-01 


2 . 5 0 OE- 0 1 




121 


122 


929 


500E-01 


2 . 50 OE - 0 1 




122 


123 


930 


500E-01 


2 . 500E-01 




123 


124 


931 


S00E-01 


2 . 500E-01 




124 


125 


932 


500E-01 


2 . 500E-01 




125 


126 


933 


500E-01 


2 . 5 0 OE- 0 1 




934 


230 


12 


500E-01 


2 . 500E-01 




935 


249 


230 


500E-01 


2 . 5 0 OE- 0 1 




936 


268 


249 


500E-01 


2 . 500E-01 




937 


287 


268 


500E-01 


2 . 500E-01 




938 


306 


287 


500E-01 


2 . 5 0 OE -01 




939 


325 


306 


500E-01 


2 . 500E-01 




940 


344 


325 


500E-01 


2 . 500E-01 




941 


363 


344 


500E-01 


2 . 500E-01 




942 


382 


363 


500E-01 


2 . 500E-01 




933 


126 


382 


500E-01 


2 . 500E-01 





LOAD CURVE DEFINITIONS 



1 2 

. OOOE+OO . OOOE+OO 
1. 000E+00 .OOOE+OO 
2 2 

.OOOE + OO 2 . OOOE-06 
1. OOOE+OO 2. 000E-06 



PRESSURE BOUNDARY CONDITION CARDS 



1 

1 

1 



17 8 

7 13 14 

13 19 20 



2 - 1 . OOOE+OO - 1 . OOOE+OO- 1 . OOOE + OO -1 .OOOE + OO 
8 - 1. 000E+0 0 - 1. OOOE + OO -1 . OOOE+OO -1 .OOOE+OO 
14-1. OOOE+OO -1. OOOE + OO -1 . 0 OOE + OO - 1 - 00 OE + 00 



FLUMAS DATA FOR BEAMSHELL MODEL 

beam, flu beam. geo stmam.dat beam.daa $ FLUNAM GEOHAM GRDNAM DAANAM 



T T F T $ 

T F T F $ 

F F T T . $ 

F F F F *' $ 

F T F F $ 

F F F F $ 

F F F F $ 

DYNA $ 

0 942 0 450 $ 

10 0 $ 

0 112 0 $ 

33 $ 

0 $ 

9.3S6E-0S 57600.0 $ 

5 $ 

1 . 1 . $ 

0 . - 1 . 0 . 0 . $ 

3 $ 

0.0 16.813 1 180 1 $ 

16.813 0.0 181 360 1 

0.0 0.0 431 520 1 

0 $ 

553 585 16 .813 16 .‘813 2 $ 

553 553 0.0 16.813 1 

553 585 16 .813 16 .813 32 



PRTGMT 


PRTTRN 


PRTAMF 


CALCAM 




EIGMAF 


TWODIM 


HAFMOD 


QUAMOD 




PCHCDS 


NASTAM 


STOMAS 


STOINV 




FRWTFL 


FRWTGE 


FRWTGR 


FRESUK 




RENUMB 


STOGMT 


ROTGEO 


ROTQUA 




PRTCOE 


STRMAS 


SPHERE 


ROTSYM 




OCTMOD 

MAINKY 


CAVFLU 


FRWTFV 


INTCAV 




NSTRC 


NSTRF 


NGEN 


NGENF 




NBRA 


NCYL 


NCAV 






NHAS 


NHAF 


NHAI 


NFUN 




NS EG 
NCIR 
RHO 


CEE 








NVEC 
CQ ( I ) 
DHALF 


CXHF 


CYHF 


CZHF 




NSRADI 

RADI 


RAD2 


J3EG 


JEND 


JINC 


NSORDER 

N1 


N2 


R1 


R2 


NSET 



AUGMAT DATA FOR 
strnam.dat beam 
F F F F 
F F F T 
F F F T 
F F F F 
F 

DYNA 

0.S0 

942 2826 3 3 
1 

0 1 450 1 
3 

0 1 86 5 

0 91 176 5 

0 400 472 6 



BEAMSHELL . MODEL 



$ STRNAM FLUNAM GEONAM PRENAM 



$ 


FRWTGE 


FRWTST 


FRWTFL 


PLNWAV 


$ 


FLUSKY 


DAAFRM 


SYMCON 


DOFTAB 


$ 


PRTGMT 


PRTTRN 


PRTSTF 


PRTAUG 


$ 


MODTRN 


STRLCL 


INTWAT 


CFADYN 


$ 


LUMPFM 








$ 


MAINKY 








$ 


DAA2 








$ 


NSTR 


NSFR 


NFRE 


NFTR 


$ 


NSETLC 








$ 


NDICOS 


JSTART 


JSTOP 


JINC 


$ 


NUMCON 








$ 


ICON 


NSTART 


NSTOP 


NINC 
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r=d3dump 

TIME INTEGRATION DATA FOR MODEL BEAMSHELL 



beam.pre beam.pos 


$ 


PRENAM 


POSNAM 






beam. rst 


$ 


RESNAM 








F F F 


$ 


BUB PUL 


REFSEC 


FLUMEN 




T F F F 


$ 


EXP WAV 


SPLINE 


VARLIN 


PACKET 


F T F F 


$ 


HYPERB 


EXPLOS 


DOUBDC 


VELIMP 


1 


$ 


NCHARG 








0.0 


$ 


HYDPRE 








0.0 -25.3125 326.4 


$ 


XC YC 


ZC 






0.0 -25.3125 16.813 


$ 


SX SY 


.SZ 
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$ 


JPHIST 








1 . 0 


$ 


PNORM 








14 . 0E-06 


$ 


DTHIST 








1 


$ 


CHGTYP 








60.0 25.8 10.0 


S 


WEIGHT 


SLANT 


CHGDEP 




1 0 


$ 


NT I NT 


NCHGAL 






0.0 2.0000E-06 


$ 


STRTIM 


DELTIM 






5000 5000 


$ 


NSAVER 


NRESET 






0 0 0 0 


$ 


LOCBEG 


LOGRES 


LOCWRT 


NSTART 


F F 


$ 


FORWRT 


STBDA2 






F 


$ 


DISPLA 
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